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Aerosol activities 
at GMAO 



□ Developing a hierarchy of global models capable of skillfully 
representing 

■ the global aerosol distribution as depicted by available in-situ and 
remotely-sensed measurements 

■ the microphysical processes needed for parameterizing cloud/ 
precipitation-aerosol feedbacks 

■ Aerosol interaction with earth-system components 

□ Developing a comprehensive aerosol data assimilation 
capability for constraining and calibrating aerosol transport models, 
including the estimation of emissions needed for driving such models 

□ Developing an aerosol forecasting capability in support of 
NASA field campaigns. 

□ Developing an aerosol observing system simulation 

capability for aiding planning of future NASA observing missions. 
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From weather to seasonal to decadal time scales 
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Hybrid EnKF/3D-Var 


NRT Forecasts 



Aerosols 


MERRAero 


EOS Instruments 


Local Displacement Ensembles 


NRT Forecasts 
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Seasonal Prediction 
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2014 NRT 
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■ 

GEOS-5 Model Configuration 


Bulk Scheme 

Dust 

Sea-salt 

Surfates 

Carbonaceous 


Finite-volume 

Cubed-Sphere 



Global, 25 km, 72 Levels, top at 0.01 hPa 
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GEOS-5 SO 2 & Sulfates 
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Non-hydrostatic 7km global meso-scale simulation 






Global, 14 km, 137 Levels, top at 0.01 hPa 
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QFED: Quick Fire 
Emission Dataset 


□ Top-down algorithm based on MODIS Fire 
Radiative Power (AQUA/TERRA) 

□ FRP Emission factors tuned by means of 
inverse calculation based on MODIS AOD data. 

□ Daily mean emissions, NRT 

□ Prescribed diurnal cycle 


JCSDA: inclusion of geo-stationary information 
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GEOS-5 Biomass Burning Aerosols ^ 



Non-hydrostatic 7km global meso-scale simulation 






Modeling Interannual Variability 
of Biomass Burning Emissions 


at 


□ BB emission anomalies respond □ 
directly to precipitation and 
surface humidity conditions 

□ The normalized Canadian Fire 
Weather Index captures the 
fammability conditions as a 
function of surface meteorology 
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Parameterization: 
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Maritime Continent Example 


Current Status of GEOS-5 Seasonal Prediction 


Application 

System Characteristics 

Decadal, IPCC 
(1960-2009) 

2 deg atmos, 1 deg ocean 
MOM4, Fortuna 2.5, specified aerosols 
(PCHEM), EnOI/SAFE(sea ice) with scout and 
MERRA forcing; original run in several streams, 
lost water - sea level problem due to precip 
correction, new run one stream/no correction 

S-l, NMME 
(1979 - present), 

1 deg atmos, >2 ocean, MOM4, Fortuna 2.5, 
specified aerosols (PCHEM), EnOI, EnKF, with 
MERRA forcing, altimeter assim. 

IESA 

(2013/14) 

% deg atmos, 1/10 deg MOM5, Latest 
Ganymed, using FAST (EnOI with time lagged 

estimates of cov) 

Current tuning 
coupled model 

2 deg atmos, 1 deg ocean, MOM5, latest 
Ganymed 

ODAS dev. 

Aquarius salinity 
analysis 

1 deg atmos, 1/2 deg ocean, MOM5, latest 
Ganymed 

Next S-l System 

>2 deg atmos. Ganymed, % (1/5) MOM5, interactive 
aerosols, new diurnal layer, updated sea ice model, 
wave model, force with MERRA2, dynamic vegetation, 
EnOI, Aquarius SSS, new synthetic salinity prior to Argo 
(2004/5), GRACE ocean mass anomalies, sea level from 
TOPEX, Poseidon, Jason 


Equatorial Pacific Temperature section 
and anomaly: AUG2011 
ENKF 



GMAO GEOS-5 Apr 2014 Forecast 



Aerosol Data Assimilation 
2014 NRT Configuration 


Focus on NASA EOS □ Simultaneous estimates of 

instruments, MODIS for now background bias (Dee and 

= Mu Siluu VJUb) = 

□ Adaptive Statistical Quality 
Control (Dee et al. 1999): 

□ State dependent (adapts to 
the error of the day) 

□ Background and Buddy 
checks based on log- 

□ Global, high resolution 2D transformed AOD 

AOD analysis 1 innova tion 

□ 3D increments by means of □ Error covariance models 

Local Displacement (Dee and da Silva 1999): 

Ensembles (LDE) □ innovation based 

□ Maximum likelihood 
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GEOS-5 Meteorological DAS 

Hybrid 3D-VAR/EnKF 


member 1 
forecast 


TT “ e n s em ;, "Analysis 


member 2 
forecast 


Generating new ensemble 

Uses background error perturbations given the latest mprn hpr 1 
covariances computed set of observations and a 

from the ensemble 

/ 

EnKF ensemble 
perturbations 
are "re-centered* 
around the high- 
res analysis 




EnKF 

member update 


first-guess ensemble 
used to estimate 
background error 
covariances 

^ * « 




forecast — 


GSI 

Hybrid Ens/Var 


Replaces the EnKF 
ensemble mean 
analysis 
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Used for GEOS-5 
Forecasts for 
Next cycle 



Deterministic 

forecast 


Previous Cycle 


Current Update Cycle 

In collaboration with NOAA NCEP, ESRL 













GEOS-5 Aerosol Assimilation 
Phase I: EnKF Only 


member 1 
forecast 


member 2 
forecast 


member 3 K : \ 
forecast 



Generating new ensemble 

Uses background error perturbations given the latest , mpm u pr * 
covariances computed set of observations and a ,tsl tutu 
from the ensemble ^ first^guess ensemb le 

EnKF ensemble 
perturbations 
are "re-centered 
around the high 
res analysis 
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s en 
esu 


background erro 
covariances 


Replaces the EnKF 
ensemble mean 
analysis 


Used for GFS 
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Aerosol Observing System 


CATS 

(ATLID) 

CALIOP 



MOD'S (NNR) 
VfIRS (NNR?) 
CEO (NNR?) 

j MISR {Bright} 
AVHRR (NNR) 
SeaWIFS (NNR?) 


LDE 


EnKF f 



Observing System 



/ AERONET 



GEOS-5 Aerosols 
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Ground 

, — — 
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Reflectances 


Forecasting 
Reanalysis only 


MGPIS 

vims 

CEO 
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Mapper 

Limb 
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POLDER 
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EnKF 


Aerosol Data Assimilation in GEOS-5 
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Global 5-day chemical 
forecasts 

• 03, aerosols, CO, 
CO2, SO2 

• Nominally 25 km 

• Driven by real-time 
biomass emissions from 
MODIS FRP (QFED) 

• Constituents 
transported on-line, 
interactively 

Since 2007 supported 
several field missions 
including TC4, ARCTAS , 
GloPac, ATTREX, 
DISCOVER-AQ, 
HS3,SEAC4RS, etc. 



Time (Hours UTC) 


Comparison of observed (top) and simulated (bottom) 
aerosol backscatter for a slight during the 2013 
SEAC4RS campaign. 


http://gmao.gsfc.nasa.gov/forecasts/ 
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NASA 


International Cooperative 
for Aerosol Prediction 
(ICAP) 


ICAP is an 

international forum 
for aerosol forecast 
centers, remote 
sensing data 
providers, and lead 
systems developers 
to share best 
practices and discuss 
pressing issues facing 
the operational 
aerosol community. 



NMMB/BSC 

Pv ■» 





ICAP Multi-Model Ensemble Dust Forecast 


Sessions et al., ACPD, 14, 14933-14998, 2014 






GEOS-5 

r\ 

NASA 

VW 

Reanalysis Activities 



Name 

Nominal 

Resolution 

Period 

Aerosol 

Data 

Available 

MERRA-i 

50 km 

1979-present 

NONE 

now 

MERRAero 

50 km 

2002-present 

MODISC5 

now 

FP for Inst. 
Teams 

50 km 

1997- 

MODISC5 

In progress 

NCA 

25 km 

2010-11 

MODISC5, 

MISR 

Now 

MERRA-2 

50 km 

1979-present 

AVHRR, 

Qi 2015 


MODISCs, 

MISR, 

AERONET 
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Observational Bias 


P\ +£ 

NASA 


ORIGINAL MODIS AOD BIAS CORRECTED AOD 





CDR r. 


AVHRR NOAA CDR AOD # 
AERONET Comparison 




Coastal and Island Stations for 2003-2009 
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r V(0 (MERRAero Analysis) 


AERONET 

MERRAero Validation 


□ 

□ 


□ 




0.99 


2.71 


2.71 
0.99 
0.36 
0 .13 
0.04 
0.01 


0.01 


0.04 0.13 0.36 

(AERONET) 


Aerosol Optical Thickness PDFJ2007-2010) 


0.56 

0.48 

0.40 

0.32 

0.24 

0.16 

0.08 

0.00 


AERONET - MERRAero 



log(r+0.01J 


rj = log(r + 0.01) 
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Regional Climatology 
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PMajj Ug/m ] PM 2 5 E^g/m 


PM 25 (Nitrates) 
Regional Climatology 
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Height (kmj Height (km) 


Vertical Structure: 
Comparison to CALIOP 




DUST 
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MERRAero - Day 



-60 -50 -40 -30 -20 -10 0 

Longitude 

CALIOP - Klinht 


10 20 -60 -50 -40 -30 -20 -10 0 

Longitude 

MERRAero - Night 



-60 -SO -40 -30 -20 -10 

Longitude 


30 -20 -10 

Longitude 


, 0,0020 


0.0015 


- 0.0010 


1 0.0005 
1 0.0000 



Evaluation of MERRAero 
Absorption using OMI 



OMI UV Aerosol Index (jul 2007) 
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Clear-Sky <% 

Aerosol Direct Radiative Effect 


AtrotOl OirKt todiMiv* tff«t (TOA) 




uo* uo*w w-m r «rt wi we- 


Source 

TOA SW DRE 

Ocean ( ) 

Atmos. 

Ocean ( id) 

Surface SW DRE 

Ocean ( ) 

MERRAero 

- 3-8 ( ) 

2.8 (6.8) 

-6.6 (-11.1) 

Other Observational 
Yu etal. (2006) 

- 5-5 ± °-2 ( ±0.7) 

3-3 (6.8) 

-8.8 ± 0.7 (-11.811.9) 

Multi-model Ensemble Yu etal. 
(2006) 

-3.4 ± 0.6 (-2.8 ± 0.6) 

( 4 - 4 ) 

-4.8 ± 0.8 (-7.2 ± 0.9) 

GEOS-5 (Free) 

- 3-4 (-2.7) 

0.5 (2.8) 

- 3-9 (- 5 - 5 ) 


DREsw = \ 

[f^sw f\w] 

1 -1 

Aerosols 

[f^sw f\w] 

1 

NoAerosols 






Assimilated Aerosol 
Annual Mean Mass 



Speciation potentially adjusted by spectral reflectances 


32 





Mass Budget 


rx l 


Annual mass budget for an aerosol specie q : 


V • (u q) = E+ P - L+ 

T 


where 

uq 

Mass flux 

E 

Emissions 

P 

Chemical production 

L 

Loss processes 

A q a 

Anaysis increments 

T 

Analysis interval (3 hours) 

(•) 

Mass weighted vertical integral 

(•) 

Time average 
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Annual Mean 
Analysis Increments 


MERftAero 2007 10 Annual Mean Dust Mass Potential y (lAUJ 





MERRAero 2007-10 Annual Mean Sulfate Mass Potential 
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O.S.S.E. 

rx 


■ Observing System 

■ Simulation 

■ Experiment 

Model-based OSSE 


A framework for numerical experimentation in which 
observables are simulated from fields generated by an 
earth system model, including a parameterized description 
of the observational error characteristics. 

Simulations are performed in support of an experimental goal. 
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Elements of an 
OSSE System 


Forecasting 

Analysis 

Retrievals 


Diagnostics 


Sensitivity 


Improvements Legacy instruments 


Calibration/Validation 


Legacy Instruments 



Applications 


Global 



t 

Regional 


f 

Meteorology 

Complexity 

Aerosols 


Chemistry 



Radiances 

Retrievals 


Physical 


Statistical 
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MODIS Cloud & Aerosol «.« 

S ^ 

Retrieval Simulator 


a) Actual RGB composite b) Simulated RGB composite 

■ PDF based sub-grid 
sampling of GEOS-5 
fields (ICA) 

■ Spatial "clumping" 

■ Radiances for 27 
MODIS channels 

■ Cloud and aerosol 
extinction, ssa, phase 
function 

■ Operational Retrievals 

■ Clouds: MOD06 
Aerosols: MOD04 

Wind et al., 2013, GMD 
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MODIS Retrieval Simulator 
Biomass Burning in Brazil 


Modis RGB Image with MYD04 retrieval overlay 


r t 



Simulated 

Retrieval 
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Modis RGB Image with GEOS-5 aerosol overlay 
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PACE Simulator 


■ GEOS-5 non- 
hydrostatic 7 km 
atmosphere with 
GOCART aerosols 

■ Coupled to GEOS-5 10 
km ocean component 
with biogeochemistry 

■ Simulation of 

Water leaving radiances 
t.o.a. reflectances 


P\ 

NASA 








Concluding Remarks 




AEROSOLS IN GEOS-5 

■ The GEOS-5 Earth Modeling 
System includes data assimilation 
of its major components 

■ Aerosols are an integral part of 
the GEOS-5 NRT and re-analysis 
systems 

■ GEOS-5 OSSE activities in 
support of new NASA observing 
missions 

Builds on NWP capabilities, extends it 
to constituents and other components 


GEOS-5 evolution 

■ Aerosol/cloud processes evolving 
from bulk to modal/2-moment 
schemes 

■ Aerosol assimilation evolving 
into a EnKF sub-system within 
the atmospheric ^D-EnVar 
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COR r. 


AVHRR NOAA CDR AOD 

MERRAero, AERONET Comparison 



MERRAero 
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AVHHR Dust Aerorot dot tea! Death (20081 



MERRAero 


AVHHR Sea-salt Aeroroi Optical Deoth (20081 
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Multiple Species 



AVHHR Sulfate Aeroroi Ontiral Deoth (20081 




PATMOS-x 

AVHRR Pathfinder Atmospheres 


PATMOS-x Dataset 

□ Version 5 Level 2B 

□ 0.1 degree sampling (not 
average) 

□ Period: 1978-2009 

□ Inter satellite calibration 
(MODIS reference) 

□ Bayesian probabilistic 
cloud detection (CALIPSO 
reference) 

■ cpd <0.5% 


- Extended 



Neural Net Retrival 

□ Ocean Predictors 

■ TOA Reflectances 

■ 630 and 860 nm 

■ TPW 

■ Ocean albedo (wind) 

■ Solar and sensor angles 

■ GEOS-5 fractional AOD 
speciation 

□ Target: 

■ AOD at 550 nm 

■ Balanced MODIS NNR 
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AVHHR Sulfate Aeroroi Optical Depth (2008} 
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Russian Super Bolide of 2013 



Gorkavyi et a!., 2013 (2013) 3:20 UTC 15 February 2013 


Chelyabinsk Meteor 
Dust Plume 




Gorkavyi eta!., 2013 (2013) As seen by OMPS Limb Profiler 





Chelyabinsk Meteor & 

Dust Plume 



Gorkavyi et al., 2013 (2013) 


As simulated by GEOS-5 


